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Within the framework of thermofield dynamics, the wavefunctions of the thermalized displaced 
number and squeezed number states are given in the coordinate representation. Furthermore, the 
time evolution of these wavefunctions is considered by introducing a thermal coordinate represen- 
tation, and we also calculate the corresponding probability densities, average values and variances 
of position coordinate, which are consistent with results in the literature. 
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I. INTRODUCTION 



Displaced number states and squeezed number states are generalization of coherent states and squeezed states 
of a harmonic oscillator, respectively 0. The coherent state is constructed by displacing the ground state of the 
harmonic oscillator p]||]j and the squeezed state by first squeezing the ground state and then further displacing it ( 
sometimes, by first displacing and then squeezing, or by only squeezing ) |Q. In these constructions, number states 
( also called Fock states in quantum field theory ) of the harmonic oscillator taking the place of the ground state 
will correspondingly produce the displaced number state and squeezed number state. Thermalizing the displaced and 
squeezed number states, one can get the thermalized displaced number state and squeezed number state, which will 
■^j- ■ be discussed in the present paper. Evidently, the thermalized coherent and squeezed states are the special cases of 
the thermalized displaced number and squeezed number states. All the above states are interesting and important in 
physics. 

As is well known, the coherent state can describe the coherent light, and its over-completeness gives rise to the 
coherent-state representation which is very useful in quantum optics, statistical physics, quantum field theory and 
particle physics, etc. |^,^|. This state mimics the motion of classical particles, and hence is also used for studying 
the Schrodinger cat states 0. For the squeezed state, not only is it a minimum uncertainty state and similar to 
the classical motion, but also the quantum fluctuations in position can be suppressed at the expense of enhanced 
fluctuations in momentum which is different from that of the coherent state. Therefore, the squeezed state has 
important technology applications in quantum computation and sensitive measurement j^j. So far the squeezed state 
has been receiving a great deal of investigations ^j,|l(| . In the same way, a displaced number state follows the motion 
q-i of a classical particle as well as keeps its shape in the course of the motion M, and the squeezed number state can 
' display the similar squeezed property to the squeezed state ( Eqs.(38) — (41) in Ref. [|J ) and hence promises hopeful 
. £h [ applications in optical spectroscopy, communications, molecular and solid state physics ( This reference dealed 
with the squeezed displaced number state ) Q. Early in 1950's, the displaced number state and the squeezed number 
state were proposed and studied with the help of the coordinate representation p2| . Since then, these states received 
a few further investigations Jl^,|lJ]. Recently, Ref. jj] reviewed these investigations, and gave the most general time- 
dependent wavefunctions and probability densities of them in the coordinate representation. Particularly, in view of 
the expremental realization of the optical and atomic squeezed ( not displaced ) states as well as the number states 
Pj (1996) [|l5|, Nieto predicted that in the not too distant future, it is hopeful to observe the displaced and squeezed 
number states 

On the other hand, no thermal noises exist nowhere, and thus the influence of the noises on the above-mentioned 
states has to be studied. Such a investigation is often realized by using density matrices and master equation. 
However, within the framework of thermofield dynamics fl6|| , a thermalizing operator acting on the states is also 
an important and useful way to introduce finite temperature effects [ |17|Jl8| . Both the density-matrices and the 
thermofield-dynamics investigations give rise to a varieties of thermal partners of the above-mentioned states, such as 
the thermalized coherent, squeezed, displaced number, and squeezed number states, the displaced thermalized state, 
squeezed thermalized state, and so on. Many properties of various thermal coherent and squeezed states have been 
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studied by constructing directly a state vector ||17|-|19| and other methods, such as characteristic function, density 
operators, Glauber's P-representation of density operator, etc. j[4 20 -2|] ( Most of Refs. po| , ^ll were concerned with 



the coherent and squeezed thermalized states ). The connections between these thermal states have been revealed in 
Ref. and Fearn and Collett also gave the physical interpretations of these states [la]. Besides, for the thermal 
coherent state, Barnett and Knight discussed the independence of the Glauber's p-representation upon the order of 
displacing and thermalizing operators [0. As for the thermalized displaced number and squeezed number states, there 
were few investigations of them, and just recently the thermalized squeezed number state ( not displaced, different 
from the state in the present paper ) was considered with its characteristic function for analysing the influence of 
thermal noise on higher-order squeezing properties of it |24| 

This paper will address the wavefunctions and position probability densities of the thermalized displaced num- 
ber and squeezed number states ( Hereafter, two of these states will also imply that the thermalized coherent and 
squeezed states are their special cases). This problem hasn't, to our knowledge, discussed in the literature, except for 
Ref. ^0|(1993) and Ref. j2^] (1965) in which the position probability density of coherent thermal state and squeezed 
thermal state ( not include number state ) was given by Glauber's R- function and/or P-represcntation ). How- 
ever, this problem is certainly interesting and meaningful. The wavefunctions of the coherent, squeezed, displaced 
number and squeezed number states contain all information about these states and hence describe completely these 
tates. Therefore, the wavefunctions of the corresponding thermalized non-classical states will give the influence of 
finite temperature on the properties described by the zero-temperature wavefunctions, and can provide, at least, a 
quantum-mechanical intuitional understanding for us. Moreover, the coordinate representation of their density op- 
erators can be obtained from the finite-temperature wavefunctions and consequently these wavefunctions can equip 
a coordinate-representation way for calculating the expectaion velues of all physical observables on the thermalized 
non-classical states, which is most usual way in quantum mechanics. Additionally, the position density probability 
can give the probability density of magnetic component of electromagnetic fields pa] Ref. p2|](1965). 

Thermofield dynamics is unique formalism of finding the wavefunctions for the thermal non-classical states. In 
this paper, within the framework of thermofield dynamics, we shall give wavefunctions of the thermalized displaced 
number and squeezed number states in terms of the position coordinate, consider their time evolution, and calculate 
the position probability densities. In order to do so, we shall first derive the wavefunction of the thermal vacuum in 
the coordinate representation, which was almost given in Ref. , and introduce a thermal coordinate representation 
in the next section. Then the wavefunctions of the thermalized displaced number and squeezed number states will 
be given in terms of the position coordinate in Section III. Section IV will address the time evolution, the position 
probability densities, the position average values and variances of these states. We will conclude this paper at the 
end. 

By the way, thermofield dynamics will be not introduced in this paper, and good expositions of them can be found 
in Ref. Besides, although this paper will discuss a harmonic oscillator with a mass and constant frequency, 

taking the mass as unit one can get the results which are usable for a one-mode electromagnetic field with the same 
frequency. 

II. THERMAL VACUUM AND THERMAL COORDINATE REPRESENTATION 

In the fixed-time Schodinger picture, for the quantum one-dimensional oscillator 

H = —p 2 + -muj 2 x 2 = (a) a + -)Tiuj , (1) 
2m 2 2 



the ground state in the coordinate representation is the wavefunction 



. muj , i . muj 9 , 
<^l°>=(^) Tea; P{-^ 2 }, (2) 



where p = —ifi-4- = —ihd x , m is the mass, u> the angular frequency, and 



-(ip + mux) , a) = (— ip + mux) (3) 



\j2mhiu ^Jlmhuj 

are the corresponding annihilation and creation operators, respectively. It is noticed that in Ref. [Q, m,uj and fi all 
are unit. In order to consider thermal effects, thermofield dynamics introduces a copy of the physical oscillator Eq.(l) 
( called the tilde oscillator ) 

H = J—p 2 + \muj 2 x 2 = (a f a + \-)Tiw , (4) 
2m 2 2 



2 



according to the tilde "conjugation": CO = C*0 |16|]. Here, C is any coefficient appeared in expressions of quan- 
tities for the physical system, O any operator, the superscript * means complex conjugation, and O represents the 
corresponding operator for the tilde system. Exploiting the physical and tilde oscillators, one can have the thermal 
vacuum [|l6| 

|O,£>=T(0)|O,O>, (5) 

where, |0, >= |0 > |0 > is the product of ground states of the physical and tilde oscillators, (3 = with fcb the 
Boltzmann constant and T the temperature, and the unitary transformation T(9) ( called thermal transformation ) 
is 

T{9) = exp{-9(f3)(aa - a^)} (6) 

with 



tanh[6(0)) 



e 



Notice that any physical operator commutes with any tilde operator. Consequently the thermal- vacuum average value 
agrees with canonical ensemble average in statistical mechanics. 

It is evident that the thermal vacuum (5) is similar to the two-mode squeezed states discussed in Ref. [|| except for a 
minus difference between the exponents in Eq.(6) here and Eq.(37) there. Although the wavefunction of the two- mode 
squeezed state was given in the coordinate representation |8j , here we still derive the position wavefunction of the 
thermal vacuum for the sake of both the completeness and the establishment of the thermal coordinate representation. 
Substituting Eq.(3) into Eq.(6), one can read 

T{9) = exp{i—(xp — xp)} (7) 

with 9 = 9{f3). From Appendix B.4 in Ref. ]26| ], the last formula can be unentangled as 

T(6) = exp{— tanh(9)xd x }exp{ln[cosh(9)](xd x — xd x )}exp{—tanh(9)xd x } . (8) 
Using the following operator properties |^7j 

e cd *f(y) = f(y + C) (9) 

and 

e Cv9 *f(y) = /(ye c ), (10) 

which are proved easily, we obtain the wavefunction of the thermal vacuum as 

™,„n .VniO^l r mill, n n., 

< x,x%P > = T(e)( — )*exp{- — (x 2 + x 2 )} 
nil 2n 

= (—)hxp{- — [(xcosh(9) - xsinhie)) 2 + (xcosh(6) - xsinh(9)) 2 }} . (11) 
irk 2n 

When (3 — * oo, < x, a;|0,/3 > is reduced to < i,a:|0,0 >. This expression Eq.(ll) can be generalized to the Gaussian 
wavefunctional approach for equilibrium field theory in thermofield dynamics p8| . 

Such an expression of the thermal vacuum wavefunction Eq.(ll) suggests the usefulness of introducing a thermal 
coordinate representation. In thermofield dynamics, for any operator of the physical or tilde oscillator Q, its thermal 
counterpart is defined as Qp = T(9)QT^ (9) |l6[] . In particular, for the fundamental canonical conjugate pairs {x,p = 
—ihd x } and {x,p = ihd x }, the corresponding thermal operators are 

x p ee T{9)xT\9) = xcosh{9) - xsinh(9), pp = T{9)pT\9) = pcosh(9) - psinh{6) (12) 

and 

xp ee T{9)xT\9) = xcosh{9) - xsinh(9), p = T{9)pT\9) = pcosh{9) - psinh(9) . (13) 
Obviously, the thermal vacuum wavefunction Eq.(ll) can be written as 
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,„ _ ,mw.i , raw, 2 -2m /has 
< x,x\O,0 >= (—)zexp{- — (xl+xl)} , (14) 

which is the same form with the wavefunction < x,a;|0,0 >. Noticing that the commutators [xp,pp] = ih, [xp,pp] — 
—ih and [Op, Op] — hold, one can set pp = —ihj^ and pp = ih-^j, and establish a representation for the thermal 
oscillator, in which any object ( operators, wavefunctions ) can be expressed in terms of xp, xp, and/or In 
this paper, we shall call it thermal coordinate representation. Evidently, this representation is reached through the 
unitary thermal transformation T{&) of the coordinate representation. When working in the representation, quantities 
will take similar forms to those in quantum mechanics, and hence it will simplify our derivation in the present paper. 

It is suitable here to mention a mathematical property and the physical sense of the thermal transformation Eq.(6) 
. It is shown easily that the action of T(9) on a function of the physical and tilde positions {a:, x} amounts to just 
the thermal coordinate xp, xp taking the place of x, x, that is, 

T(9)f(x,x) = f(xp,xp). (15) 

Thermal transformation is also called thermalizing operator Jl8| . It describes the effect of a thermal reservoir in which 
a quantum harmonic oscillator immerses. From Eq.(5), we can say loosely that a thermalizing operator heats the 
ground state of a zero-temperature harmonic oscillator into a thermal vacuum with a finite temperature. In quantum 
optics, the thermalizing operator describes the action of a source which excites one-mode electromagnetic field from 
its ground state to a chaotic state ( thermalized radiation ). Thus, in order to consider thermal noise, it is enough to 
perform the action of the thermalizing operator on the non-classical states mentioned in the last section. Next, we 
shall address them. 



III. THERMALIZED DISPLACED NUMBER AND SQUEEZED NUMBER STATE IN THE 

COORDINATE REPRESENTATION 

Because both coherent and squeezed states are constructed with the displacing operator and squeezing operator 
acting on the ground state, there are three different states with squeezed effect : squeezed state ( only the squeezing 
operator acting on the ground state ), displaced squeezed state, and squeezed displaced state, which are all usually 
called squeezed state in the literature. In this paper, the terminology "squeezed state" means only the displaced 
squeezed state, for which the action of the displacing operator follows that of the squeezing operator. So does the 
squeezed number state. However, when introducing a finite temperature effect, one still faces more choices about 
the orders among displacing, squeezing and thermalizing. A different order will lead to a different thermal non- 
classical state |18|j . Nevertheless, if using thermal creation and annihilation operators to work, i.e., doing as done 
in Ref. ]l7| , ^9[ , one can escape the order problems with thermalizing operator. In this section, we shall introduce a 
finite temperature effect into the displaced number state and squeezed number state by using the thermal creation and 
annihilation operators with the vacuum |0, (3 > |l7| , |l9| ] and then give their expressions in the coordinate representation. 
This construction is utterly to thermalize the displaced number and squeezed number states, namely, it gives the 
thermalized displaced number and squeezed number states, as one shall be seen later. 

The thermal annihilation and creation operators with the thermal vacuum Eq.(5) are [ fl6| 

ap = T(0)aT\G),ap = T{Q)a)T\9) (16) 

and 

hp = T{6)hT\6),h p = T{6)tfT\9) . (17) 

One can easily check that ap\0,[3 >= 0, ap\O,0 >= and [0^,0^] = [ap, at] = 1. With the aids of thermal creation 
operators ap and at, one can construct normalized thermal number states 

\n, m, (3 >= -jt=a\ "4 m |0, > (18) 

with the closure relation 

\n, m, >< 0, m, n\ = 1 . (19) 

n,m 
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The so-called displaced number state \a, n > of the oscillator Eq.(l) is defined in Fock space as \a, n >= D(a)\n > 
and can have the following form in the coordinate representation [l[ 

ttilu i 1 . 
<a;a,n> =( — — 1 4 . — = ea^H— zaia^l 



r / 2ft 2 2muj jmuj r- . . . 

' ex P\ — 7rr\ x -\ a V +i\ —^a2x}H n [\ ——x~ V2ai 20 

2n V mw V n \ n 

with a = (ai + ia.2) any complex number, H n [- ■ •] the Hermite polynomials and the displacing operator 

D(a) = e aaf - a ' a . (21) 

In this definition, when the number state \n > is replaced by the ground state |0 >, the state \a, n > is reduced to the 
usual coherent state \a >. Evidently, the state \a,n > is constructed just with the displacing operator D(a) acting 
on the number state \n >. Similarly, one can define the following state |a, n, (3 > 

\a,n,/3 >= D^(a)^(a)|n,n,/3 > (22) 

so as to introduce a finite temperature effect into the displaced number state. Here, the thermal displacing operators 
D p{a) and Dp{cy) are 

Dp(a) = exp{aaL — a*ap} (23) 



Dpipi) = exp{aap — a*ap} , (24) 

respectively, which are generalization of the displacing operator D(a). Note that a = a* in the present paper ( Of 
course, one can take a as another parameter independent of a ). When n = the state \a, n, (3 > is just Eq.(ll) with 
7 = a in Ref. [^9)( the first paper ) and Eq.(3.1) with ip = a in Ref. jL7). Employing the definitions (5), (16), (17), 
(18), we obtain 

\a, n, (3 >= T(6)\a, n>\a,n> . (25) 

This equation indicates that the state \a,n,(3 > is just the thermalized displaced number state. In the last equation, 
\a, n > is the tilde displaced number state and can be obtained from Eq.(20) according to the tilde rules. When n = 0, 
the state |a,0,/3 > is the thermalized coherent state, being similar to Eq.(3.3) in Ref. p7| . Employing Eqs.(15),(12) 
and (13), one can have 



~ 1 n ,mw.i 1 r muj,. / 2h ,0 / 2k 2 

<x,x\a,n,(3> = — 2 ^-^ exp {- X0 _ W Ql ) + _ W ai ) 

7m 2™rt! 2n \ ttilu V 



+«y — — a 2 (a; / 3 - a; /3 )}ii rl [W — x/3 - V2aiJii„[W — xp - V2«iJ (26) 



mw 1 1 tow _ . / 2?L 2 



+ (xcosh(9) — xsinh(d) — \l a i) 2 } + i\l —T~ a 2( c °sh(8) + sinh(9))(x — i)} 

" ttilu V n, 



i\J-J[-( xcosh ( ) - xsinh(d)) - V2ai]H n [^^-(xcosh{e) - xsinh{6)) - V2«i] . (27) 



The r.h.s of Eq.(26) is the wavefunction of the thermalized displaced number state in the thermal coordinate repre- 
sentation, and Eq.(27) is just the wavefunction in the coordinate representation. When (3 — > 00, Eq.(27) is reduced 
to a product of i-function and i-function, each factor resembling Eq.(15) in Ref. Q. 

Now we are at the position to discuss the thermalized squeezed number state. The squeezed number state \a, z,n > 
of the oscillator Eq.(l) is constructed by using the squeezing operator S(z) 0] 



S(z) — exp{~ — (z* aa — za) a))} (28) 
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and reads 

\a, z, n >= D(a)S(z)\n > . (29) 

Here, z is any complex constant. When n = 0, |a, z, > is the usual squeezed state. From Ref. |5J], the wavefunction 
of \a, z, n > in the coordinate representation is ( here in terms of our notations ) 

i ,mw i (a/Ts)" r . , 

<x\a,z,n> = — — 4 — = . exp\—iaia,2\ 
1 v 7rft ; V-^i2™n! 



mw 2h 2 2muj jmuj / 2?i 

•ea;p{-— ^2 x - \/ ai + — — a 2 a;}ii„U/ — (J^J {x-\ ai) (30) 

Ah V muj \ n y n V mu> 

with z = zi + iz2 = re 1 ^ 1 S = cosh(r) + z\sinh{r)/r and k = Z2sinh(r) / (2rS) . In Eq.(30), for the convenience of 
comparison later, we adopted the notations T's in Ref. H, that is, 



T 1 =S{1 + i2k), T 2 = ; — -%2k, 



1 

S 2 {l + i2n) 



r* = hr%r> ^ = s(i + ^ . (31) 

In analogy with the definition of the squeezed number state, we introduce the thermal squeezing operator 

Sp(z) = exp{-^(z* apap - za^a^)} , Sp(z) = exp{~(z*apap - za^a^)} (32) 

and define the following state \a, z,n, f3 > 

\a, z, n,/3 >= Dp(a)Dp(a)Sp(z)Sp(z)\n, n,0 > . (33) 

to introduce a finite temperature effect. Note that z = z* in the present paper ( Of course, one can take z as another 
parameter independent of z ). It is easily shown that 

\a, z, n, (3 >— T(9)\a, z,n > \a, z, n > (34) 

with \a, z,n > the tilde version of \a, z, n >. Evidently, the last equation indicates that the state \a, z,n,/3> is just 
the thermalized displaced number state. When n = 0, \a, z, 0, /3 > is just the thermalizd squeezed state, being similar 
to Eq.(ll) in Ref. flg]( the second paper ). Employing Eqs.(15),(30),(12) and (13), we obtain 



mui^i | 71 , mw 2Ti 2 



< x,x\a,z,n,(3> = { — j^7—-exp{- — [T 2 {x p - \ a x ) 

irn LFi \2 n n\ 2n V raw 



/ 27i n2 i 2muj _ 
-Fzixp - \l — ax) \ + «y— ^— a 2 (xfi - xp)\ 



•ffn ^-T ^ ^ - V a lPn \/ -Z-(^4)~ (5/3 " \/ — «l) 35 

1 a V tow V ft V mw 



+^2 (xcosh(9) — xsinh{9) — \l Oi\) 2 ] 

muu 



+i\ — r — a 2 (cosh(9) + sinh(9)){x — x)} 



•H n \il 1 (xcosh(9) — xsinh(9) — \j — — ai) 

n V tow 



■Iln^l^iF^ixcoshie) - xsinh(6) - \j ^-ax)] . (36) 



J7JW 



G 



The expression Eq.(36) is just the wavefunction of the thermalized squeezed number state in the coordinate represen- 
tation. When (3 — > oo, the x-part of Eq.(36) is consistent with Eq.(20) in Ref. 0. 

In this section, we have constructed the thermalized displaced number and squeezed number states with the thermal 
creation and annihilation operators, and given their wavefunctions. These states are physically meaningful. For a 
displaced thermalized squeezd state, Fearn and Collett gave a physical interpretation that it corresponds to the output 
from a linear photon amplifier whose input is a squeezed state if the amplifier's added noise is regarded as thermal 
photons Thus, according to this interpretation, it is not difficult to give physical interpretations for the states 
in the present paper : thermalized displaced number and squeezed number states. For instance, the thermalized 
displaced number state should correspond, at least theoretically, to the output from a thermal source who excites the 
previous output from a linear photon amplifier with a number state being the input, if no thermal noises company 
the amplifier and the thermal source can excite an electromagnetic field from its ground state to a thermal chaotic 
state. Of course, strictly speaking, it is impossible to have no thermal noises, and so one should consider a type of 
completely thermalizied state in which both the displacing and the squeezing are companicd by thermal noises. We 
shall discuss the more practical situations in a separate paper. Next, we shall discuss the time evolution of the two 
thermalized non-classical states. 

IV. TIME-EVOLUTION OF THE THERMALIZED NON-CLASSICAL STATES 

In this section, we first consider the time-evolution of the thermalized displaced number and squeezed number states 
and then calculate position probability densities of them. 

In thermofield dynamics, Hamiltonian H of the combined system of the physical and tilde oscillators is jl6| 

H = H — H = (a) a — aJa)hu> = (agdp — aXap)hu> . (37) 

Hence the time-evolution operator of the combined system is Jl6[ | 

U(t) = exp{——Ht} = exp{~iujta ii l3 ap}exp{iujtaj 3 ap)} (38) 
with t the time. From Eqs.(3), (16) and (17), we have 

ap = 1 (ipp + muxp) , at = 1 (-ipp + mujxp) (39) 
v2mfiw y2mnuj 

and 

~ap = 1 {-ipp + mujxp) , at = 1 (ipp + muxp) . (40) 
\2mnuj y/2mnu! 

Thus, in the thermal coordinate representation, the time-evolution operator can be uncntanglcd as 

, , 1 r muj 2 . .. <9 r h d 2 
U (t) = — -exp{— i— — tan{uit)xp\exp{— log(cos(ujt))xp— — \exp{i- tan{u!t)—-i^\ 



Tfioj d Ti 

■exp{i-—tan(uot)x 2 p\exp{—log(cos(ujt))xp-^—\exp{—i- tanitut)— y }. (41) 

2n ' oxp 2mu) ax^ 



The operator U (t) acting on a wavefunction will yield the time-evolution of the wavefunction. With the help of Eqs. (9) 
and (10) and the following operator property j27| 

exp{Cd 2 y }f(y) = -j= J exp{-^fL}f(w)dw , (42) 

which can be shown by using the identity 

i />oo 

v2tt J-oc 

we can perform the action of U (t) on various thermalized non-classical states here. 
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Leting U(t) Eq.(41) act on the wavefunction Eq.(14), one can find that the wavefunction Eq.(14) is invariant and 
hence the wavefunction of the thermal vacuum is independent of time, i.e., < x,x\0,/3, t >= U(t) < x,x\0,/3 >=< 
x,x\0,(3 >. This is understandable because the average value of any physical observable on the thermal vacuum is 
equal to its ensemble average value, which does not vary with time. 

U(t) Eq.(41) acting on Eq.(26) gives the time-dependent wavefunction of the thermalized displaced number state 
< x, x\a, n, /3, t >= U(t) < x,x\a,n,(3 > as 

muj l 1 a 2 a* 2 . , . 2l 

< x,x\a,n,p,t > = ( — — ) 2 - — rexp\ (— H — — )cos{ujt) — 2a, \ 



■ex P {-^[{ X0 - ^ — -jf + (a* - y — ^) 2 ]} 

■H n [iS—^-xp — V2(aicos(cut) + a2sin(u>t))]H n [^J -jr-Xp — V2(aicos(ujt) + a>2sin(L>t))] (43) 

with A — cos(tut) + isin{ujt). In finishing the relevant integral with Hcrmitc polynomial, we used the formula 7.374 
(8) in Ref. |29j . Eq.(43) indicates that the wavefunction of the thermalized displaced number state is dependent upon 
time. 

Similarly, U(t) Eq.(41) acting on Eq.(35) yields the time-dependent wavefunction of the thermalized squeezed 
number state < x, x\a, z, n, (3, t >= U(t) < x, x\a, z,n,j3 > as 

< x, x\a, z, n, p, t > = ( — —J 



irh' \F x \2 n n\\B\ 

, T2Cos(u>t)a\ + 2T2sin{uj)oi\oi2 + isin(uJt)oa , 
■exp{ } 

, TXcos(uit\oi\ + 2TX sinUS)a.\a.2 — isinhot)oa , 
■ ex P{ -jp 1 



muj T2cos(ujt) + isinkot) 2 Imuj F^oti + ict2 

■ exp{ ~^h b 'e + 2 V 2» b X0} 



mu> J^cosiut) — isin(u>t) 2 lmu> T^ai — ia,2 

' exp{ ~^n W* X0 + 2 V ~2h w* X0} 



/ raw is— 1/ / 2fi 

U/ -r ^ B (acfl - y {cos{u!t)ai + sz7i(wi)a 2 )) 

V ft V wo; 



■^nU/^^l-Bl)" 1 ^ " \/— (coa(wi)ai + «n(wi)a 2 ))] ■ (44) 

V ft V TflUJ 

with i? = cos(ujt) + i!F2sin(ujt) . Like the thermalized displaced number state, the wavefunction of the thermalized 
squeezed number state is also time-dependent. Substituting Eqs.(12) and (13) into Eqs.(43) and (44), one can obtain 
the time-dependent wavefunctions of the thermalized displaced number and squeezed number states in the coordinate 
representation. When /3 — * oo, the x-part of Eq.(44) is consistent with Eq.(45) in Ref. M, except for the lack of an 
imaginary exponent, which is cancelled by the relevant exponent from the tilde system (In Eq.(46) of Ref. [jjj there 
should be the symbol "=" between the parentheses and the fraction Y Additionally, for the case of (3 — > oo, the x-part 
of Eq.(43) is consistent with the zero-z resultant of Eq.(45) in Ref. Q] ( except for the lack of an imaginary exponent 
)■ 

Eqs.(43) and (44) indicate that the wavefunctions of both the thermalized displaced number state and the ther- 
malized squeezed number state are dependent upon time, which is different from the thermal vacuum wavefunction. 
This point is because both the thermalized displaced number state and the thermalized squeezed number state are 
not eigenstates of the Hamiltonian H, while the thermal vacuum is an eigenstate of the Hamiltonian H with zero 
eigenvalue. 

Now we can calculate the position probability densities. The probability density is the modulus square of the 
position-coordinate wavefunction with the tilde coordinate integrated. At the first, we consider the thermal vacuum. 
The density is easy to calcuate and the result is 



/oo />oo 
< t,(3,0\x,x >< x, x\0,/3, t > dx = I < (3, 0\x, x >< x, x\0, (3 > 
-oo J —oo 



dx 



raw ,i./3huj. . raw ,.3huj. . , , 

—tanh^(—^-)exp{-—tanh(—^-)x 2 }. (45) 
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This is just the familiar result in statistical mechanics. 

Secondly, we calculate the position probability density of the thermalized displaced number state 

/oo 
< t, (3, n, a\x, x >< x, x\a, n, (3,t > dx . (46) 
-oo 

Substituting Eqs.(12),(13) and (43) into Eq.(46) and reducing it, we have 

Pc(x,t)^^-(—— f / ex P {-( ai x + a 2 ) 2 - {b 1 x + b 2 ) 2 }{H n [a 1 x + a2} f(H n [-b 1 x-b 2 }) 2 dx (47) 

where, 



ai = 



Imoj FrrwJ . r-. . . . 

W -j—cosh(9), a,2 = — d -j—xsinh(9) — v2(aicos{ujt) + a 2 sin(ujt)) 



bi = iJ -y—sinh(8), b 2 = —^—xcosh(6) + y/2{a,\cos{ut) + a 2 sin(uit)) . 
With the help of the formula on page 225 of Ref. fl3(| 

min{m,n} 

H m [x]H n [x}= 1 r r\C r m C r n H m+n _ 2r [x] (48) 

r=0 

with C T n and being combinations, Eq.(47) can be written as 

Pc(x,t) = ^(^) 2 £ V +k 3 \k\{CiCl? \ exp{-y 2 - (ay + b) 2 }H 2n ^ 3 [y]H 2n - 2k [ay + b]dy (49) 
nwi z n. jk=Q j_ x 

with a = ^- and b = —aa,2 + b 2 - Using repeatedly the formula 

e- x2 H n [x] = -±{e~^ H n ^[x]} 
and then the formula 7.374(8) in Ref. p9[, one can obtain 



mui 1 2 / 7r r b , 

Pc(x,t) = ——(——) J exp{ T —\ 

nnai 2 n n\ y a + 1 a + 1 

• V +k mCiC^) 2 a 2 ^(a 2 + iy+"-^H2 {2n - 3 - k) [^—} (50) 



•ea;p{ — [y — —tanh\ )x — \Z2(a.\cos(wt) + ot2Siri(uit))y 1 + secfe( - )j } 
• £ 2 2 ^^)- 2 " J U!(^^) 2 e^{i(^fc)^/3}(co S M^)) J+ ' £ - 2 ' i 



j,fe=0 



r mljJ sfiftu. rr, , . . / f3flUJ 

■-ti2(2n-j-k)[y —r-tanh(—^-)x - V 2(a.iCOs(Lut) + a 2 sm(ujt)) y 1 + sech( — — )\ (51) 

As for the thermalized squeezed number state, the calculation of the position probability density is completely 
similar to that of the thermalized displaced number state. Finishing calculations similar to the above, one can have 
the position probability density of the thermalized squeezed number state Eq. (44) 
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p s (x, t) = I (|a, z, n, /3, £ >)* | a, z, n, /3, £ > (52) 

J — oo 

muJ 7 4 \B\ t 1 2 / i (.F3.F3)" _ r fo 2 



" 7 

• 2*+*jU!(C£C*)V<»-JV + xy+^ng [ ° ] (53) 

^ x/a 2 + LF 4 |5| 

mw / 1 ^ 2 1 j. l 1 / /3fiw ^ r 1 



tan) ^FT^r tan/l2 (-^-) ea; p{ 



tt/i v 2«nK Ti\B\ y 2 ' l v T 2 B*B 



, mUJ rr, , s . / \\ A, 7//3^. l2l 

•[y -j—tanh{ - )x — v 2(a:icos(u;t) + a2S«n(wi)j y 1 + sech{ ^ )j ) 

• £ 2 2 ^+ fc )- 2 ™ J U!(^^) 2 e^{i( J -fc)^/3}(c OS M^)) J+fc - 2 " 

1 /mcj Bhuj . r-. , . , ,. / 

7777 (y -^—tanh( )x — v2(aicos{ujt) + a.2Sin(u)t))\ 1 + sech(— — ))\ . (54) 
i> \ n 2 V 2 



i,k=o 



'"2(2n-.j-k) I 

When /3 — * oo, the existence of the factor a 2 (" _J ) enforces the summation index j have a unique value n because 
a = 0. Thus, employing Eq.(48), one has 

E 2 i+fe j!fc!(C^)V^V + ly+fe-^g [ ] = 2"n!(g n [ ]) 2 , 

3 ^Z Va 2 + 1^4 1 -B| Va 2 + LF4I-BI 

and hence Eq.(54) with 3 — > oo can give Eq.(50) in Ref. jjj. Meanwhile, the probability density Eq.(51) with 8 — > oo 
can also lead to Eq.(50) with z = in Ref. Q. 

For a given time, for example, t = 0, one can get the position probability densities p c (x) and p s (x) without 
considering the time evolution from Eqs.(51) and (54). Making a contrast between p c (x) and p c (x,t) as well as 
p s {x) and Ps(x, t), one can find that just the displacement and squeeze parameters in the expressions of the densities 
experience changes with the evolution of time. That is, for the thermalized displaced number state, the real part of 
the displacement parameter a becomes ot\cos(u)t) + ia2sin{tot) , which is similar to that in Ref. p2]| ( 1965 ), and for 
the thermalized squeezed number state, besides the same change of «i, the parameter T± becomes T^\B\. 

The expressions of both p c (x,t) and p s (x,t) are complicated, but one can easily prove that each of them is nor- 
malized, because only the term with j = k = n is not zero when p c (x,t) or p s (x,t) is integrated with respect to x. 
Furthermore, one can calculate the average value of the position coordinate x on the thermalized squeezed number 
state as 



< x >= / xp s (x,t) = \/cMh( ^ W )\ ^^(aicos(ujt)+a2sin(ujt)) . (55) 
J-oo V 4 V tow 

Evidently, < x > is independent of n and the squeeze parameter z, which is similar to the result of the squeezed 
number state Q. Eq.(55) indicates that the average value of x on any thermalized squeezed number state follows 
the motion of a classical harmonic oscillator, and the amplitude of the oscillation increases with the increase of the 
temperature. When 8 — > oo, Eq.(55) is consistent with Eq.(36) in Ref. and contrasting to Eq.(36) in Ref. [Q, 

Eq.(55) has just an additional temperature factor J coth(^^-) . 

Besides, one also easily obtain the variance of i on a thermalized squeezed number state as 

(A n x) 2 >^ ^ (x- < x >) 2 p s (x, t) = coth(^)(2n + . (56) 

When 8 — > oo, (A n x) 2 is consistent with Eq.(38) in Ref. |l]]. When n — 0, (Aqx) 2 is consistent with Eq.(15a) in 
Ref. [fL9|( the second paper ). A comparison of Eq.(56) with Eq.(38) in Ref. jy] tells us that (A n x) 2 just has an 
additional temperature factor coth{@^). 

Finally, we give two examples of p s (x,t) to end this section. For n = 0, we have 



, ^ / 1 r (x- <x >) 2 , r _. 

Ps{x,t) = , —- — ^exp{ — - — ^— . (57) 

W 2n{A x) 2 2(A x) 2 
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Using the relation between the thermalized squeezed state and squeezed thermalized state in Ref. E^)(1991), one can 
find that the last equation with t — is consistent with Eq.(6.6a) in Ref. |^(|(1993). For n = 1, we have 



/ 2 (a- < a; >) 2 f l 2 fffio; (a:- < x >) 2 3 2 3 1 (x- <£>)% 

^ s(x ' t} = y^W 6a:P{ ~ 2(A x)^ { 2 Secfe ( — )[( 2(A ^ - 2 } 2 ] + 2(A x)^ 1 " (58) 

From the last equation, we see that by introducing a finite temperature effect, x-polynomial factor in the expression 
of p s (x, t) is not just the quadratic term of [x— < x >) as Eq.(52) in Ref. pj, but have additional quadruplicate term 
of (x— < x >). The appearance of the quadruplicate term is understandable because thermalizing a nonclassical state 
amounts to doubling the freedom number of the systems within the framework of thermofield dynamics Jig ] . 

From Eq.(55) to the last equation, we give some explicit results only about the thermalized squeezed number state. 
As for the thermalized displaced number state, taking z = in Eqs.(55) — (58), one can obtain the corresponding 
results about them, which are also consistent with those in the literature. 



V. CONCLUSION 



In this paper, we have given the wavefunctions of the thermalized displaced number and squeezed number states 
in the coordinate representation. Furthermore, with the help of the thermal coordinate representation, we obtain the 
time-dependent expressions of these wavefunctions. Although the thermal vacuum wavefunction is time-independent, 
but either the thermalized displaced number state or the thermalized squeezed number state varies with time. We 
also give the probability densities, average values and variances of the position coordinate on these stats. Each of 
the wavefunctions, the time-dependent wavefunctions, the probability densities, average values or variances here is 
consistent with that in the literature when the temperature tends to zero. Setting n = in the expressions of this 
paper, one can obtain results of the usual thermalized coherent and squeezed states. Additionally, setting t = 0, one 
can obtain the probability densities without considering the time evolution of these states. 

In the thermal coordinate representation, the forms of the wavefunctions Eqs.(14),(26),(35), (43) and (44) resemble 
their own zero-temperature limits. Of course, this resemblance does not exist in the coordinate representation at all, 
and the probabilty densities are different from those at the zero-temperature case. From section IV, one has seen 
that the thermal coordinate representation simplified greatly the calculations there. Perhaps, the thermal coordinate 
representation would simplify other calculations related to the thermal non-classical states. 

Finally, although the thermalized displaced number and squeezed number states are discussed, the above results 
of them can give easily the corresponding ones of other similar states, such as displaced thermalized number state, 
squeezed thermalized number state, etc., with a simple parameter transformation JT^] (1991) J2Q] (1993). Additionally, 
no matter how complicate the expressions (51) and (54) are, it is not difficult to compute them numerically for a given 
number n. We believe that once the displaced number state and squeezed number state is prepared in laboratories 
some day, the results in the present paper will be useful. 
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